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The finite size effects of the thermal conductivity k have been studied in the phonon space. It is 
found that only a few phonon modes are selected to take part in the thermal transport when the 
size L of the system is decreased. The amount of the selected phonon modes is proportional to the 
L. In this way, k decreases with the decreasing of L. Such mechanism for the size effect of k found 
in this work is beyond the Phonon-Boundary scattering. The exponent a of the power law k ~ L a 
has been fitted, showing that the exponent is not universal. 

PACS numbers: 44.10.+i, 66.70.-f, 65.80.-g, 63.20.-e, 63.22.-m. 


I. INTRODUCTION 

With the development of nano technologies, electronic 
devices can be fabricated with scales down to mesoscopic 
size. In such scale, some of the physical properties 
of mesoscopic system change dramatically with the 
system size. For example, while thermal conductivity 
k for dielectric or semiconductor materials is a size 
independent in macroscopic system, it depends on 
the system size drastically in mesoscopic scale 00 - 
The size effects on the n have been confirmed by the 
experiments B0. Especially the experiment basing on 
the transient thermal grating (TTG) shows that the k 
increases with increasing system size and is converged 
to a finite value for the infinite system size 0. There 
exist two contradictions between such experimental 
results and the known theories. One is that calculations 
based on the harmonic interaction of atoms give rise 
to a divergent k even if the disorder and anharmonic 
interaction in the model are involved [Tol-[T^. This is 
the so-called divergence problem in thermal transport. 
The second contradiction is the decreasing of n with 
the system size. Some authors have introduced a soft 
cut-off l c for the mean free path of phonons 
They claimed that the phonons with the mean free path 
l > l c contribute little to the thermal transport. But 
how l c appear is still unknown. In this paper, we give a 
simple model which can resolve above two contradictions. 

Generally, the thermal transport for one dimensional 
system is simplified as the coupling between oscillators 
ordered in chain. For a chain with a finite length, 
the oscillator number is noted by L and the lattice 
parameter is a. The heat baths are attached to the two 
ends of the oscillator chain as a whole model dMl. it 
has been rigorously proved that if the coupling between 
the oscillators is harmonic, the thermal conductivity n 
diverges with the increasing of the size La EH. Even 
by invoking the disorder and anharmonic interaction 
in the model to enhance the thermal resistance, k is 
still not converged for the infinite size (TH, Hi: which 
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is the well-known divergence problem of k 00. The 
divergence problem has been revealed by some other 
models in the real space, such as the hydrodynamics 
model [16], and the model basing on the Mode-Coupling 
theory 17]. According to those models, the relation 
between the k and L follows a power law by n ~ L a with 
a > 0. The exponent a varies from model to model 00. 

In the phonon space, the divergence problem is 
naturally removed by the phonon scattering. The 
phonon scattering enhances the thermal resistance and 
converges the k to be a finite value kb ~ cv 2 t for the 
infinite size [HI, [HMEU] . The phonon scattering has been 
implemented in the relaxation time r. Here, c is the 
heat capacity of the system and v is the averaged group 
velocity of phonon. The mechanisms for the phonon 
scattering are various, such as by Umklapp scattering 
EH, Iffl. Molecular Dynamics(MD) simulation has 
also been used to study the mechanism [26i-[29|. In 
order to simplify the calculation, the approximation 
of single-mode relaxation time Tk has been applied to 
identify the contribution of each phonon at a wave 
vector k in the Brillouin Zone(BZ). Then the kb is 
modified to be kb ~ f ckVi. Tkdk with the integral over 
the whole BZ [H [l8l422i j. Unfortunately, no entire 
model covering the whole BZ is available for the Tk 
yet. Some authors have proposed that Tk oc 1 [uj 2 for 
phonon scattering at high frequency uik [H, E0 El- 
Such relation can not be extended to the low frequency 
range. Or, it is unphysical that Tk goes to infinity 
when ujk is approaching to the T point of BZ, which 
also leads to an divergent kb■ In order to simplify the 
phonon scattering for low frequency phonon, we make a 
assumption to modify the model by Tk oc l/(w^ + C ) for 
the extension. Here, C a positive constant, dominating 
the Tk when tok goes to zero. The introduction of 
C into the Tk is reasonable since the mean free path 
vt of phonon has been measured to be a finite value [32| . 

The size effect on the n can be obtained in the term 
of Phonon-Boundary(PB) scattering 12714291 . [33:]. 

The PB scattering means that phonons are scattered 
by the boundaries with the relation of r ~ L, which 
decreases the k for the finite system size j28|. But this 
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is in fact incorrect in TTG case since the TTG is formed 
within a Si thin film and no any interfaces exist for the 
PB scattering (§]. In the TTG case, the authors have 
introduced a soft cut-off l c for the mean free path of 
phonons to show the size effect on the n @,[ 0 . The l c 
means that the phonons with the mean free path l > l c 
contribute little to the thermal transport. But why does 
the TTG system have the l c ? In this work, we will show 
that only a few phonon modes are selected to take part 
in the thermal transport instead of the total phonons in 
the BZ. Such mechanism is beyond the PB scattering. 


II. THEORY 

We consider the TTG as an uniformly ordered chain 
in one dimension with an infinite length. The oscillator 
number in one grating period is the L. The interaction 
between the oscillators is harmonic and the phonon 
scattering has been implemented in the relaxation time 
Tfc. The phonon dispersion relation of the TTG can be 
obtained by using the lattice dynamics theory and the 
wave vector k is obtained by applying the Born-Von 
Karmen Boundary Condition(BKBC) on the chain 0- 
The total Hamiltonian H of the system is the summation 
of the Hamiltonian of each oscillator by H = )T/ m H rn ■ 
Here, H m is the Hamiltonian of the m th oscillator, 
which is H m = p 2 m /(2M) + (1/2) J2 n x m$mnX n - In this 
expression, x m is the displacement of the rn th oscillator 
deviating from its equilibrium position and p m is the 
momentum. <f> mra is the force constant for the interaction 
between the m th and the n th oscillators. And M is the 
mass of one oscillator. 

According to the Heisenberg’s equation of motion, 
the time derivative of H m is H m = J2 n {®™nX m Pn - 
d ) nm^nPm]/(2Af). On the right hand side of the equation, 
the first term means the work done by all the other oscil¬ 
lators on the m th oscillator while the second term is for a 
reversed process. We introduce the creation operator a/" 
and annihilation operator a*, for phonon at k , and rewrite 
the total Hamiltonion as H = J2k( a k a k + 1/2 )hx k - The 
operators of x m and p m can be expressed as functions of 
Ofc and a/", and then the time derivative of H rn reads 

H m = Y,[ e - i(k+r)RmT kAt) + Wk, r (t)\ + h.c. 

k,r 

(i) 

Here, T'k, r and W k , r are time-dependent operators, as 
functions of a k ( r ) and a k ^ in the phonon space. R m is 

the equilibrium position of the m th oscillator. 

Since the energy of each oscillator is real and periodical 
in the chain by H m+ L = H m , it can be obtained from the 
eq. m that k or r takes only the values of 

A =-£—, h = 0,±1,±2,---,±—. (2) 


That means only the phonon at A are selected to take 
part in the time evolving of H m and are responsible 
for the decaying of the TTG. All the other phonon 
modes not at A in the BZ act as a constant thermal 
background, and are eliminated after the time deriva¬ 
tive. For the selected phonon, the nearest wave vector 
close to the T point is k c = ±2ir/(La) with h = ±1, 
indicating that phonon in the range of 0 < |fc| < \k c \ 
contributes little to the thermal transport even though 
with their large mean free paths. Thus, k c acts as the 
origin for the_soft cut-off l c proposed in the models for 
the TTG [§, 0. The amount of the selected phonon 
modes A equals L. That means fewer phonon modes are 
selected for the thermal transport for smaller L by mod¬ 
ulating the thermal grating. In this way, the size effect 
of the k can be clarified. This issue will be specified later. 

Before we calculate the n for the TTG, we turn to an¬ 
other type of the oscillator chains with a finite length 
of the total chain, comparing to the TTG with a finite 
length of one period but an infinite length of the total 
chain. Such finite-size cliain(FSC) can be found like the 
nanotubes [§J. We will get the similar result of the TTG 
for the FSC. The FSC is a uniformly ordered chain in one 
dimension with two open ends, and L now represents the 
total number of the oscillators of the chain. For conve¬ 
nience, we index the ordered oscillators from 0 to L — 1. 
The Hamiltonian of the chain reads 

Hf=Y, MX ™ Xm + Xm^nXnQmOn- ( 3 ) 

m m,n 

Here, 0 m („) is the geometric structure factor(GSF) of the 
chain, and defined as 0 m („) = 1 for 0 < m(n) < L—l and 
0 m („) = 0 for other cases. By using the GSF, the FSC 
is extended to be infinitely long without any changing to 
the total Hamiltonian. The motivation for such extension 
is that the lattice dynamics theory and the BKBC can 
only be applicable on the chain with an infinite length. In 
this way, we can get the basic phonon dispersion relation 
as a start point for the study of FSC. We use N to de¬ 
note the total number of oscillators of the extended chain 
with TV — > oo. We rewrite the Hamiltonian by using the 
Fourier transformation x m = ( l/y/NM ) e lkRm Qk■ 
The Qk is the normal coordinate of phonon at k for the 
extended chain. The Hamiltonian then reads 

H f = \Y, QkQ*Zk- r + A q Z q _ k Z r _ q QlQ r . (4) 

k,r k,r,q 

Here, Z k is the GSF defined in the phonon space, which 
is obtained by the Fourier transformation of 0 m . The Z k 
has an expression of Z k = ^^ e ika_i) • The basic phonon 
dispersion relation u> q ~ q of the extended chain is 
obtained by the Fourier transformation of the force con¬ 
stant 4> mn . It can be found that the normal coordinates 
Q are coupled to each other due to the finite size effect. 
It is also found that when L is approaching to TV, the 
GSF Z k goes to S k ,o and eq.@ recovers the diagonalized 
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Hamiltonian H f = \YlkQkQ* k + hHk^lQlQk f° r a 
real infinitely long chain. 


We introduce a new set of coordinates S in the phonon 
space by the definition of E p = QrZ r - P ■ The motion 
equation of the S p then reads 


d 2 Z p a r /a 2 e iio ( fe -P) - 1„ „ „ 

dt 2 + 2 jt J_„ /a e iQ ( fc -rf - 1 




in integral form. The rigorous numerical solution to the 
integral equation is not applicable due to the highly os¬ 
cillated integral kernel. It is expected that the density 
distribution of the phonon eigenmodes as the solutions 
are narrow peaks located at some wave vectors. We ap¬ 
proximate the density peaks as 6 functions and denote 
the locations of the eigenmodes by p or q in the follow¬ 
ing. Thus, the density of states (DOS) of p mode phonon 
is noted by 2Trr] p 5(k —p). The DOS is used to transform 
the eq. ([31) to be in the form of summation, reading 


e iLa(q—p) _ i 

n p ^p =a2_^U 2 eia{q _ p) _ 1 2-qVq ( 6 ) 

9 

with Q p the frequency of the p mode and E p ~ 

The phonon eignmodes of the solution should be de¬ 
coupled to each other, meaning q — p = 2nh/(La) must 
be held to reduce the eq.© to be = aLr] p uj 2 E p . 

Here, h represents integers. Then, one obtains that 
the locations of the phonon eigenmodes p(q) of the 
FSC are at A = 2tt h/(La) with h the integers in the 
range of [—L/2,L/2). This result is the same as that 
of TTG, only different in meaning of L. When L is 
approaching to N for the infinitely long chain, the 
equation of motion for the p mode should recover 
S p = — WpSp. That means the eigenfrequency fl p must 
be equal to oj p and r/ p equals 1 /(La). Physically, 2nr] p 
is the degeneracy of the p mode per unit length of the 
FSC. It is interesting that applying the BKBC directly 
on the FSC by using the plane waves as the lattice waves 
can get the same result of A. However, it should be 
emphasized that the direct application of the BKBC is 
not rigorous for the study, since the DOS of phonon eign¬ 
modes are not exactly the 5 functions according to eq. ©. 


Fig.1(a) shows the phonon dispersion relation of 
the FSC by the discrete data points while the solid 
continuous line is the relation for the infinitely long 
chain. For the calculation in this paper, the dispersion 
relation of c <j(k) = wo| sin(fc/2)| is used as an example 
with wo = 4 x 10 12 Hz of Silicon applied [Hj. And 
the lattice parameter is set to be a = 1 for simplicity. 
Comparing to the solid line, phonon modes not at A are 
absent for the FSC. Especially, the phonon in the range 
of — 2n/(La) < k < 2ir/(La) have no contribution to 
the thermal transport for the FSC. The absence of the 
phonon close to the T point for the thermal transport 
has been confirmed by the TTG experiment [fj. 
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FIG. 1: (a) Phonon Dispersion Relations(PDR). Solid 
and continuous line is the PDR for the infinitely long 
chain while the discrete data points are for the chain 
with a finite size. In the plot, L = 20 is used. (b)The 
renormalized thermal conductivity k/kb as a function 
of the system size L with various C. In the plot, C have 
been renormalized by Wq. (c) In (k/kb) as a function of 
ln(L), showing the power law of k oc L a for small L. 

The exponent a can be fitted from the curves. 


III. RESULTS AND DISCUSSION 

We have studied the phonon modes of two types of 
the mesoscopic structures. One type of the structures 
is the TTG and the other type is the FSC. The two 
types of the structures have the same result that only a 
few phonons located at A can take part in the thermal 
transport. We call these phonons by selected phonons. 
The total amount of the modes of the selected phonons 
equals L. The period length of TTG or the total length 
of the FSC, both of which have been denoted by La 
in this paper, can be changed to bring three effects on 
the phonon behaviors. Suppose the L is decreased as 
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TABLE I: The exponent a of the power law k ~ L a 
fitted from the curves in figl.(c) with various C. The 
error for the fitting has been indicated. Results show 
that the exponent a is not universal. 


c 

0.01 

0.001 

0.0001 

a 

1.167 

1.324 

1.34 

error (±) 

0.0486 

0.03247 

0.03017 


understand the size effect on k. In order to study the 
exponent a for the power law k ~ L a , we have plotted 
the \h(k/kb) as a function of ln(L) in fig.1(c). The 
exponents a are fitted from the curves in the range of 
6 < L < 20, and the results are listed in table 1. It 
shows that a is not universal, and dependents on the 
detail phonon scattering from material to material. 


an example to specify the three effects. The first effect 
is that the mode amount L of the selected phonons 
then is reduced in the BZ. The second effect is that 
the 77 = 1 /{La) is enhanced for each selected phonon 
mode. Actually, these two effects are canceled for the 
thermal transport due to the produce of the mode 
amount and the degeneracy as a constant of Lij = 1/a. 
The third effect is that the selected phonon mode close 
to the T point are shifted to higher frequencies by 
w = wosin( 7 r/(La)). According to the Bose-Einstein 
distribution f k (T) = ( e h ^/(k B T) _ 1 )-i ; higher fre _ 
quency state has fewer averaged phonon amount at a 
fixed temperature. Thus, the third effect leads to fewer 
phonons taking part into the thermal transport, which 
is the key mechanism to the size effect on k. 


IV. CONCLUSION 

As a conclusion, we have studied the finite size effect 
on the thermal conductivity k in the phonon space. 
The k decreases with the decreasing of the size of 
the structures. The selected phonon modes have been 
obtained to interpret the size effect, which is beyond the 
Phonon-Boundary scattering. The exponents a of the 
power law k ~ L a have been fitted, showing that the 
exponent is not universal. 


Appendix A 


Now we jump to the calculation of k. The general 
expression of the n reads 

k = 2tt^2 V p c p v 2 p t p . (7) 

p—A 

Here, 2irr] p has been understood as the degeneracy 
of the selected phonon mode. For the FSC, we have 
obtained y p = 1 /(La). And for the TTG, r] p is a product 
of two factors by r] p = 771 - 772 - Here, 771 = 1 /(Na) is 
contributed from the infinite length Na of the TTG, 
while 772 = N/L is the amount of the grating periods 
as copies of the thermal transport in the real space. 
Thus, r] p = 1 /(La) is the same for the TTG. And, c p 
is the heat capacity for the p phonon mode while 77 p c p 
is the specific heat due to the rj inversely proportional 
to La. The c p has the expression of hjj p dfk{T)/dT. 
The phonon group velocity v p has the definition of 
v p = du}(k)/dk\k= P . The relaxation time t p is still 
within the single-mode approximation, and has been 
modified to be oc l/(w(i + C) as discussed above. 
In the calculation, the summation is over the discrete 
modes A as obtained in eq. ©■ By setting L = N, the 
r] p then acts as the infinitesimal in the BZ, transforming 
the eq. 0 in integral form to recover kb ~ f CkV^T^dk. 

Fig.1(b) is the plot of k/kb as a function of L with 
various C. In the plot, C has been renormalized by 
Wg for clarity. The plot indicates the size effects on 
the k, which matches the behaviors of n measured by 
experiments SI- The k converges to kb for L —> 00 . 
The mechanism of the selected phonons is beyond 
the Phonon-Boundary scattering, and can be used to 


In this appendix, we complete the derivation of eq.(l) 
in the main text. The Hamiltonian of the m th oscillator 
reads 


Hm = + \ X] x m§m S X n - (Al) 

n 

According to Heisenberg equation, one obtains the fol¬ 
lowing equation after some algebra 


Hm — \H : Hm\ — ^ ^ 2 ^n.mXnPm] • 

n 


We define the Fourier transformations 

(A2) 

^ VNM^ e 

n Qk 7 

Q k — 

i p -ikRr 

V N ^ 

m 

n %m ■> 

k 

' Pk ; 

Pk = \ 

1 1 \ ^ r ikRrr 

1 MN ^ 

m 

‘ Pm 1 

(A3) 


with Q and P the coordinates and momenta in the 
phonon space. R m is the equilibrium position of the m th 
oscillator. Since the x and p are real, one gets Ql = Q-k 
and P£ = P-k- By using the create and annihilation 
operators 


Ofc = 



a 


+ 

k 



(A4) 
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we have 


Hu, 


a relation between the force constant $ and the phonon 
frequency oj k by 


Qk = ]J^( a k+ a t k ); Pk = i\j '—^—{ a k~ a - k ), (A5) 
Substituting ea. (IA5D into the ea. (IA3D . one obtains 


1 - M " 2 
/ J c ^mn — 


Vn 


VN“ k ’ 


^ mn — Y^ Afu; fe e 


k(R m -R n ) 


(B3) 


-i(kR m -ui k t) + 


2 N A'Iuj k 


a k (®)]j We define Zfc by 


(A6) 




AR m k 


N 


^ gikLa _ ^ 

TV e ifea _ I ■ 


Pm — 


_ __ (B4) 

V i\l [e- <(Mm -"*■ t] a+ (0) _ e dfc« m t) afc (°)]. 

fc V 2TV Actually, is the geometric structure factor in the 

(A7) phonon space. The Hamiltonian then reads 


Then the motion equation ea. (IA2l) reads 

-i {k+ r)R mTkr ^ + e^-r)R mWkr ^ + ^ 

(AS) 


Hm = y~l[e~ 

k.r 


with 


T k ,r(t) = 


W k , r (t) = 


Hi R.o. 


— [w? - l vl]e i< ' Urt+ulkt) n+(n\n+ 


4TV/TVV w fc 

Hi fijj. 


«+(0)a+(0); 

(A9) 


4A/TV V Wfe 


? - W 2] e i(uvt- Wfe t) 


a fc (0)o+(0). 

(A10) 


H f=J2 ^T7 eiRm{r+k) QkQr 


m,k,r 


2 TV 


^ e i(rfln+fca ro )Q fc Q f 


2NM 


(B5) 


m.n.k.r 


— Zr-kQ k Qr + Y, WqZ q _ k Z r _ q Q k Q r 


k,r k,r,q 

This is the eq.® in the main text. The Lagrange func¬ 
tion C is 

£=\\Y Zr-kQlQr - Y U$Z g - k Zr- q QlQr ]■ (B6) 


k,r 


k,r,q 


Appendix B 


The conjugate momentum P q of the Q* in the phonon 
space can be obtained as 


In this appendix, we complete the derivations of eq. ® 
and eq.® in the main text. The ordered oscillators of the 
finite chain are indexed from 0 to A — 1. The Hamiltonian 
of the FSC in the main text is the eq. ©, repeated here 
by 

\ Mx m Xm 1 \ - „ „ /ni'i 

■Hf — / v ^ ®m 2 / v (Bl) 

m m,n 

By using the Fourier transformation x m = 
^=== P kRm Qk, one rewrites the Hamiltonian 
as 


0, 


Hf=Y ^ iRm(r+k) QkQr 


+ Y ^§^e i{rRn+kRm) QkQr- 


(B2) 


m,n,k,r 


2NM 


nr* i 

P q = „ x = n YXQkZ k - q + Q k Z- q - k ) = Y^ Q k Z k — q . 

d Qq 2 Y V 

(B7) 

And then the equation of motion reads 

QjZj- _ 

Pq = = ~Y. U} kQrZ k - q Z r -k- (B8) 

Define a new set of coordinates by 

: YQ k Z k - p . (B9) 

k 

Combining ea (IB7l) . ea (IB8ll and ea (lB9l) . we obtain the 
equation for the coupled phonons: 

+ ’Y', u} 'kZk—pZ‘k — 0- (BIO) 


Here, TV is the total number of the oscillators of the ex¬ 
tended chain with TV — > oo, which has been specified in 
the main text. For an infinitely long chain, there exists 


Since the wave vector k and p are in the Brillouin Zone of 
the infinitely long chain, the ea. dBlOl) can be transformed 
in integral form leading to eq.® in the main text. 
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